High-energy behavior of the nuclear symmetry potential in asymmetric nuclear matter 
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Using the relativistic impulse approximation with empirical NN scattering amplitude and the nu- 
clear scalar and vector densities from the relativistic mean-field theory, we evaluate the Dirac optical 
potential for neutrons and protons in asymmetric nuclear matter. From the resulting Schrodinger- 
equivalent potential, the high energy behavior of the nuclear symmetry potential is studied. We find 
that the symmetry potential at fixed baryon density is essentially constant once the nucleon kinetic 
energy is greater than about 500 MeV. Moreover, for such high energy nucleon, the symmetry po- 
tential is slightly negative below a baryon density of about p — 0.22 fm -3 and then increases almost 
linearly to positive values at high densities. Our results thus provide an important constraint on 
the energy and density dependence of nuclear symmetry potential in asymmetric nuclear matter. 

PACS numbers: 21.65.+I, 21.30.Fe, 24.10.Jv 



I. INTRODUCTION 

Recently, there is a renewed interest in the isovector 
part of the nucleon mean-field potential, i.e., the nu- 
clear symmetrypotential, in isospin asymmetric nuclear 

matteriaiaaaiaaa&aiaiii.iiiiiiiiiiiiiia. 

Knowledge on the symmetry potential is important for 
understanding not only the structure of radioactive nu- 
clei and the reaction dynamics induced by rare isotopes, 
but also many critical issues in astrophysics. Besides 
depending on the nuclear density, the symmetry po- 
tential also depends on the momentum or energy of 
a nucleon. Various microscopic and phenomenological 
models, such as the relativistic Dirac-Bruechner-Hartree- 
Fock (DBHF) 0, H H [H [H and the non-relativistic 
Bruechner-Hartree-Fock (BHF) pl^| approach, the rel- 
ativistic mean-field theory based on nucleon-meson in- 
teractions |12j , and the non-relativistic mean-field theory 
based on Skyrme-like interactions 0, , have been used 
to study the symmetry potential. However, the results 
predicted by these models vary widely. In particular, 
while most models predict a decreasing symmetry poten- 
tial with increasing nucleon momentum albeit at different 
rates, a few nuclear effective interactions used in some of 
the models lead to the opposite conclusion. Thus, any 
constraint on the momentum and density dependence of 
the symmetry potential is very useful. 

In the optical model based on the Dirac phenomenol- 
ogy, elastic nucleon-nucleus scattering is described by the 
Dirac equation for the motion of a nucleon in a relativis- 
tic potential. For spherical nuclei, good agreements with 
experimental data were obtained in the relativistic ap- 
proach with a scalar potential (nucleon scalar self-energy) 
and the zeroth component of a vector potential (nucleon 
vector self-energy), while the standard non-relativistic 
optical model using the Schrodinger equation failed to 
describe simultaneously all experimental observables |17| . 



Motivated by the success of the Dirac phenomenology, a 
microscopic relativistic model based on the impulse ap- 
proximation (RIA) O E El was developed, and 
it was able to fit very well the data from p+ 40 Ca and 
p+ 208 Pb elastic scattering at nucleon energies of both 
500 and 800 MeV. A nice feature of the RIA is that it 
permits very little phenomenological freedom in deriving 
the Dirac optical potential in nuclear matter. The basic 
ingredients in this method are the free invariant nucleon- 
nucleon (NN) scattering amplitude and the nuclear scalar 
and vector densities in nuclear matter. This is in contrast 
to the relativistic DBHF approach, where different ap- 
proximation schemes and methods have been introduced 
for determining the Lorentz and isovector structure of 
the nucleon self-energy 0, H, IS EiJ E3| • 

In the present work, we evaluate the Dirac optical 
potential for neutrons and protons in asymmetric nu- 
clear matter based on the relativistic impulse approx- 
imation by using the NN scattering amplitude deter- 
mined by McNeil, Ray, and Wallace [22], which has 
been shown to be valid for nucleons with kinetic energy 
greater than about 300 MeV. The high energy behav- 
ior of the nuclear symmetry potential from the result- 
ing Schrodinger-equivalent potential is then investigated 
without adjustable parameters. We find that the nuclear 
symmetry potential at fixed density becomes almost con- 
stant for nucleon kinetic energy greater than about 500 
MeV. For such high energy nucleons, the nuclear sym- 
metry potential is further found to be weakly attractive 
below a nuclear density of about p = 0.22 fm~ 3 but to 
become increasingly repulsive when the nuclear density 
increases. 

The paper is organized as follows. In Section [HI we 
briefly review the relativistic impulse approximation for 
nuclear optical potential and the relativistic mean-field 
model for nuclear scalar and vector densities. Results 
on the relativistic nuclear optical potential and the nu- 
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clear symmetry potential in asymmetric nuclear matter 
are presented in Section 11111 A short summary is then 
given m Section EH 



II. NUCLEAR OPTICAL POTENTIAL 

A. Relativistic impulse approximation 

Many theoretical studies have suggested that nucleon- 
nucleus scattering at sufficient high energy can be viewed 
as the projectile nucleon being scattered from each of the 
nucleons in the target nucleus. One thus can describe the 
process by using the NN scattering amplitude and the 
ground state nuclear density distribution of the target 
nucleus. For the Lorentz-invariant NN scattering ampli- 
tude, it can be written as 

F = F s +F v ^ l2t ,+F T a^a 2 ^+Fp 1 fll+F Al h^l 1 2^ 

(1) 

in terms of the scalar Fs, vector Fy, tensor Ft, pseu- 
doscalar Fp, and axial vector Fa amplitudes. In the 
above, subscripts 1 and 2 distinguish Dirac operators in 
the spinor space of the two scattering nucleons and 7's are 
gamma matrices. The five complex amplitudes Fs, Fy, 
Ft, Fp, and Fa depend on the squared momentum trans- 
fer q 2 and the invariant energy of the scattering nucleon 
pair, and they were determined in Ref. |22T | directly from 
the NN phase shifts that were used to parameterize the 
NN scattering data. For a spin-saturated nucleus, only 
the scalar (Fs) and the zeroth component of the vector 
(-FV7172) amplitudes dominate the contribution to the 
optical potential. In the relativistic impulse approxima- 
tion, the optical potential in momentum space is thus ob- 
tained by multiplying each of these two amplitudes with 
corresponding momentum-space nuclear scalar ps (q) and 
vector /5y(q) densities, i.e., 

^opt(q) = ~^ h [F S (q)ps(<l) + loF V (q)p V (q)}, (2) 

where puh and M are, respectively, the laboratory mo- 
mentum and mass of the incident nucleon. The optical 
potential in coordinator space is then given by the Fourier 
transformation of J7 op t (q) , similar to the "tp" approxima- 
tion used in non- relativistic impulse approximation |22| . 

The q-dependence of the relativistic NN amplitude is 
important for calculating observables of a nucleon scat- 
tering off finite nuclei within the Dirac phenomenology. 
In the present work, we are interested in the strength of 
the Dirac optical potential of nucleons in infinite nuclear 
matter. Since the scalar and vector densities in coordi- 
nate space are constant in infinite nuclear matter, they 
are delta functions in momentum space, i.e., ~ 5^ 3 ^(q). 
In this case, only the forward NN scattering amplitudes, 
i.e., F S o = F s (q = 0) and F V o = F v {q = 0), contribute 
to the Fourier transform of Eq. @, and the nuclear 
coordinate-space optical potential takes the simple form 



m 

U op t = 4 ^J lab [Fsops + loFvoPv] , (3) 

where ps and py are, respectively, the spatial scalar and 
vector densities of an infinite nuclear matter. 

The optical potential J7 op t is a 4 x 4 matrix in the 
Dirac spinor space of the projectile nucleon and includes 
following scalar U s ot and vector [/q * pieces: 

U opt = U s ot + lo U t Q ot . (4) 

Since U s ot and C/q * are generally complex, they can be 
expressed in terms of their real and imaginary parts, i.e., 

U s ot = U s + iW s , Ul ot = U + iW , (5) 

B. Nuclear scalar densities 

To evaluate the RIA optical potential for neutrons and 
protons, we need the values for the isospin-dependent Fso 
and Fyo, which can be found in Ref. j22|, as well as the 
scalar and vector densities for neutrons and protons. For 
the latter, we determine them using the relativistic mean- 
field (RMF) theory |2^| with a Lagrangian density that 
includes the nucleon field ip, the isoscalar-scalar meson 
field a, the isoscalar-vector meson field u>, the isovector- 
vector meson field p, and the isovector-scalar meson field 
6, i.e., 

C(i>, a,,p,5)=i> [ 7fl (id" - g u <J>) - (M - g a a)\ i> 
+\(d^a - mla 2 ) - \uj^ v + \mlu^ 

-^b a M(g a af - ^c IJ (g a a) i + ^(g 2 , uj^f 

+ \(d^5 - m 2 5 S 2 ) + \m 2 p p^ - \p^ v 

+\(g 2 P Pn-p' l )(^sglv 2 + Avgl^) 

-gpPv^Til) + gsS^JTip , (6) 

where the antisymmetric field tensors lo^ and p^ are 
given by uj^ = - d^ijj v and p^ = d v p p - d t _ l p v , 

respectively, and the symbols used in Eq.JfJJ have their 
usual meanings. The above Lagrangian density is general 
and allows us to use many presently popular parameter 
sets. In the present work, we use three typical parame- 
ter sets, namely, the very successful NL3 model |24|, the 
Z271v model, which was used to study the neutron skin of 
heavy nuclei and the properties of neutron stars [2^| , and 
the HA model which includes the isovector-scalar meson 
field 5 and fits successfully some results calculated with 
more microscopic DBHF approach |26|. 

In Fig. 2] we show the neutron and proton scalar den- 
sities ps as functions of the baryon density pB (vector 
density in the static infinite nuclear matter) in nuclear 
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FIG. 2: (Color online) Energy dependence of real and imag- 
FIG. 1: (Color online) Neutron and proton scalar densities inary partg of the gcalar and yector potentials for neutrons 

as functions of baryon density in nuclear matter with isospin and protong fa nudear matter with igogpin asymmetry a = 

asymmetry a = and 0.5 for the parameter sets NL3, Z271v, and Q g for the parameter set HA . 

and HA. 



matter with isospin asymmetry a — and 0.5 for param- 
eter sets NL3, Z271v, and HA. The isospin asymmetry 
is defined as a = (p n — p p )/pb with pb = p n + p P and 
p n and p p denoting the neutron and proton densities, re- 
spectively. It is seen that the neutron scalar density is 
larger than that of the proton at a fixed baryon density 
in neutron-rich nuclear matter. While results for differ- 
ent parameter sets are almost the same at lower baryon 
densities, they become different when p-g, > 0.25 fm~ 3 
with Z271v giving a larger and NL3 a smaller p§ than 
that from the parameter set HA. The real and imaginary 
parts of the scalar potential at higher baryon densities 
thus depend on the interactions used in evaluating the 
nucleon scalar density and have, therefore, large uncer- 
tainties. 



III. RESULTS 

A. Relativistic nuclear optical potential 

With neutron and proton scalar densities obtained 
from the RMF theory for the parameter set HA, we have 
studied both the energy and density dependence of the 
real and imaginary parts of the scalar and vector po- 
tentials for neutrons and protons in nuclear matter with 
isospin asymmetry a = and 0.5. In Fig. |2 the resulting 
energy dependence is shown for the three nucleon densi- 
ties p B = 0.08 fm" 3 (panel (a)), 0.16 fm" 3 (panel (b)), 
and 0.24 fm -3 (panel (c)). For all densities, the optical 
potential shows a strong energy dependence below 300 
MeV, where it is known that the influences due to ambi- 
guities in the relativistic form of the NN interaction, the 
exchange contribution, and the medium modification due 



to Pauli blocking are important. The low energy behav- 
ior of the optical potential can in principle be studied in 
the generalized relativistic impulse approximation based 
on the relativistic meson-exchange model of nuclear force 
and using the complete set of Lorentz-invariant NN am- 
plitude [23,123,123,133,133. Since many theoretical studies 
have shown that data on elastic nucleon-nucleus scatter- 
ing can be reproduced by using above optical potential 
when the nucleon kinetic energy is greater than about 
300 MeV and that this optical potential also agrees very 
well with that extracted from phenomenological analysis 
of the nucleon-nucleus scattering data |18lll9ll20ll32l |. we 
thus focus in present work on the higher energy behavior 
of the isospin-dependent optical potential. As shown in 
Fig. for all three densities considered here, there is 
a systematic difference or isospin splitting in the optical 
potentials for protons and neutrons in asymmetric nu- 
clear matter. Specifically, the neutron exhibits a stronger 
real but weaker imaginary scalar and vector potentials in 
neutron-rich nuclear matter. Furthermore, both the pro- 
ton and neutron optical potentials become stronger with 
increasing density. 

The density dependence of the real and imaginary 
parts of the scalar and vector potentials for neutrons and 
protons in nuclear matter with isospin asymmetry a = 
and 0.5 obtained with the parameter set HA is shown 
more explicitly in Fig. [21 for the three nucleon kinetic 
energies of E kin = 600 MeV (panel (a)), 800 MeV (panel 
(b), and 1000 MeV (panel (c). An isospin splitting of the 
nucleon optical potential in asymmetric nuclear matter 
is again clearly seen. 
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FIG. 3: (Color online) Density dependence of the real and 
imaginary parts of the scalar and vector potentials for neu- 
trons and protons in nuclear matter with isospin asymmetry 
a = and 0.5 for the parameter set HA. 



B. Nuclear symmetry potential 

From the Dirac optical potential, a "Schrodinger- 
equivalent potential "(SEP) of the following form is usu- 
ally introduced [13, 13 : 



EP 



t/£ ot 
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tot2 



7" /"tot 



where e is the nucleon kinetic energy. Using the SEP 
in the Schrodinger equation gives the same bound-state 
energy eigenvalues and elastic phase shifts as the solu- 
tion of the upper component of the Dirac spinor in the 
Dirac equation using corresponding Dirac optical poten- 
tial. The real part of the SEP is then given by 



Re((7 ; 



SEP ) 



Us + U Q - 



1 

2M 



(Ul- 



U 



(8) 

Above equation corresponds to the nuclear mean-field po- 
tential in non-relativistic models |35| and allows us 
to obtain the following nuclear symmetry potential, i.e., 
the so-called Lane potential 36] : 



"SVm 



Re((7sEp)n - Re(c7sEp)p 
2a 



(9) 



where Re(<7sEp)n and Re((7sEp)p represent, respectively, 
the real part of the SEP for the neutron and proton. 

In Fig. we show the energy dependence of the nu- 
clear symmetry potential for the parameter sets NL3, 
Z271v, and HA at fixed baryon densities of p-g, = 0.08 
far 3 (panel (a)), 0.16 fur 3 (panel (b)), and 0.24 fur 3 
(panel (c)). It is seen that all three parameter sets give 
similar nuclear symmetry potential for a nucleon at ki- 
netic energy higher than about 300 MeV, i.e., it first 
decreases with nucleon kinetic energy and then becomes 
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FIG. 4: (Color online) Energy dependence of the nuclear sym- 
metry potential using the parameter sets NL3, Z271v, and HA 
as well as from the phenomenological interaction MDI with 
x — —1, 0, and 1 at fixed baryon densities of ps = 0.08 far 3 
(a), 0.16 fm -3 (b), and 0.24 fm -3 (c). 



essentially constant when the nucleon kinetic energy is 
above about 500 MeV. Specifically, the nuclear symme- 
try potential starts from about MeV at lower density of 
Pb = 0.08 fm -3 (about half of nuclear saturated density), 
4.8 MeV at normal nuclear matter density (pe = 0.16 
fnr 3 ), and 12 MeV at higher density of p B = 0.24 fm -3 
(about 1.5 time nuclear saturated density) and then sat- 
urates to about -3.8 ± 0.5 MeV, -1.8 ± 1.7 MeV, and 
5.3 ± 3.8 MeV, respectively, when the nucleon kinetic en- 
ergy is greater than about 500 MeV. The uncertainties 
in the saturated values simply reflect the variation in the 
energy dependence of the symmetry potential at high en- 
ergies. 



For comparison, we also show in Fig. 0]results from the 
phenomenological parametrization of the momentum- 
dependent nuclear mean- field potential, i.e., MDI inter- 
action with x — —1, 0, and 1. In the MDI interaction, 
the single nucleon potential in asymmetric nuclear mat- 
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ter with isospin asymmetry a is expressed by 0, 0, 0, IT(i| 
C/(p,a,p,T,r) 



-95.98-*^^ 



-120.57+ £ 



o" + 1 / Po 
2B \ p T 



(T + 1J po 



+ 



Po 
Po 



a-l 



+ 1 Po 



-ap T 



d 3 p' 



/r(r,p') 



d 3 p' 



1 + (p - p') 2 /A 2 
/r'(r,p') 



1 + (P-P') 2 /A 2 ' 



(10) 



In the above r = 1/2 (—1/2) for neutrons (protons) and 
t 7^ t'; cr = 4/3; and / T (r, p) is the phase space distribu- 
tion function at coordinate r and momentum p. The 
parameters -B, CV,r, C T . T > and A were determined by 
fitting the momentum-dependence of U(p, a, p, r, r) to 
that predicted by the Gogny Hartree-Fock and/or the 
Brueckner-Hartree-Fock (BHF) calculations the sat- 
uration properties of symmetric nuclear matter and the 
symmetry energy of 31.6 MeV at normal nuclear matter 
density po = 0.16 fm~ 3 The incompressibility K of 
symmetric nuclear matter at po is set to be 211 MeV. The 
different x values in the MDI interaction are introduced 
to vary the density dependence of the nuclear symmetry 
energy while keeping other properties of the nuclear EOS 
fixed [T^j. We note that the energy dependence of the 
symmetry potential from the MDI interaction is consis- 
tent with the empirical Lane potential at normal nuclear 
matter density and low nucleon energies and has been 
used in the transport model for studying isospin effects 
in intermediate-energy heavy ion collisions induced by 
neutron-rich nuclei [a, 0, ITol ] . 

It is seen from Fig. 21 that results from RIA at lower 
density of p — 0.08 fm -3 are comparable to those from 
the MDI interaction with x = 0, while at higher baryon 
density of pe = 0.24 fm -3 they are comparable to those 
from the MDI interaction with x = — 1. At normal nu- 
clear matter density, the MDI interaction, which gives 
same results for different x values by construction, is seen 
to lead to smaller nuclear symmetry potential at high nu- 
cleon kinetic energies compared with present results from 
the RIA based on empirical NN scattering amplitude and 
the nuclear scalar density from the relativistic mean-field 
theory. We note that our results agree surprisingly well 
with those of DBHF by Fuchs et ai. HE3|. 

For the density dependence of the nuclear symmetry 
potential using the parameter sets NL3, Z271v, and HA 
at a fixed high nucleon kinetic energy of 800 MeV, it 
is shown in Fig. together with corresponding results 
from the MDI interaction with x = —1, 0, and 1. It is 
clearly seen that the nuclear symmetry potential from all 
parameter sets NL3, Z271v, and HA changes from nega- 
tive to positive values at a fixed baryon density of about 
Pb = 0.22 fm -3 and then increases almost linearly with 
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FIG. 5: (Color online) Density dependence of the nuclear 
symmetry potential using the parameter sets NL3, Z271v, and 
HA as well as from the MDI interaction with x — — 1, 0, and 
1 at a fixed nucleon kinetic energy of 800 MeV. 



baryon density. Furthermore, the nuclear symmetry po- 
tential depends not much on the choice of the parameter 
sets NL3, Z271v, and HA. At such high nucleon kinetic 
energy, the nuclear symmetry potential from the MDI in- 
teraction with x = reproduces nicely the results from 
the RIA when ;$ 0.1 fm -3 as in its energy dependence 
at low densities shown in Fig. 0] The two differ strongly, 
however, at high densities. The MDI interaction with 
both x = — 1 and 1, on the other hand, show very differ- 
ent density dependence from present RIA results. 



IV. SUMMARY 

We have evaluated the Dirac optical potential for neu- 
trons and protons in asymmetric nuclear matter based on 
the relativistic impulse approximation with empirical NN 
scattering amplitude and the scalar and vector densities 
from the relativistic mean-field theory. We find that the 
nuclear symmetry potential derived from the resulting 
Schrodinger-equivalent potential is not very sensitive to 
the parameter sets used in the relativistic mean-field cal- 
culation, particularly at low densities and high nucleon 
energies, although the latter give very different nuclear 
scalar densities at high baryon densities in both symmet- 
ric and asymmetric nuclear matters. Furthermore, the 
nuclear symmetry potential at a fixed density becomes al- 
most constant when the nucleon kinetic energy is greater 
than about 500 MeV. For such high energy nucleon, our 
study shows that the density dependence of its nuclear 
symmetry potential is weakly attractive at low densities 
but become increasingly repulsive as nuclear density in- 
creases. Results presented in present study thus provide 
important constraints on the high energy behavior of the 
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nuclear symmetry potential in asymmetric nuclear mat- 
ter, which is an important input to the isospin-dependent 
transport model |j El in studying heavy-ion collisions 
induced by radioactive nuclei at intermediate and high 
energies. They are also useful in future studies that ex- 
tend the Lorentz-covariant transport model [38l to 
include explicitly the isospin degrees of freedom. 
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